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Abstract 


Typical  synthetic  aperture  radar  (SAR)  images  are  two-dimensional,  providing  range 
and  azimuth  information,  but  furnish  few  details  with  respect  to  elevation.  One 
approach  to  extend  SAR  to  three-dimensional  imaging  is  to  replace  the  single  antenna 
element  by  an  array  of  antenna  elements.  The  simplest  such  implementation  of  this 
would  be  a  linear  array  oriented  vertically.  This  report  contains  two  sections,  the  first 
outlining  data  and  image  processing  for  a  three-dimensional  application,  and  the 
second  part  details  finite  beamwidth  processing  for  the  time-domain  backprojection 
method.  Standard  backprojection  implicitly  assumes  infinite  beamwidth,  whereby  the 
entire  target  area  is  illuminated  and  returns  are  collected  from  all  points  in  the  imaging 
grid.  This  decreases  the  efficiency  of  the  algorithm  by  summing  over  unnecessary  grid 
points.  Finite  beamwidth  processing  is  designed  to  be  more  physically  realistic  by 
accounting  for  the  actual  antenna  pattern.  Initial  results  show  that  depending  on  the 
width  of  the  effective  beam,  the  number  of  processing  loops  can  be  reduced  by  up  to 
50%,  while  still  maintaining  good  image  quality  in  terms  of  the  reconstructed  target 
response  (pulse  width).  Potential  improvements  in  3D  SAR  processing  are  even 
greater. 


Resume 


Les  images  typiques  de  radars  a  synthese  d’ouverture  (SAR)  sont  bidimensionnelles  et 
foumissent  des  donnees  de  distance  et  d’azimut,  mais  peu  de  details  sur  le  site.  Un 
moyen  d’ajouter  une  troisieme  dimension  a  l’imagerie  SAR  consiste  a  remplacer 
l’element  d’antenne  unique  par  un  reseau  d’elements  d’antennes.  La  configuration  la 
plus  simple  serait  un  reseau  lineaire  oriente  verticalement.  Le  rapport  se  divise  en  deux 
parties;  la  premiere  porte  sur  le  traitement  des  donnees  et  des  images  dans  une 
application  tridimensionnelles,  et  la  seconde  examine  en  detail  le  traitement  de 
largeurs  de  faisceau  finies  selon  la  methode  de  retroprojection  dans  le  domaine 
temporel.  La  methode  de  retroprojection  courante  suppose  l’usage  d’antennes 
omnidirectives,  au  moyen  desquelles  l’ensemble  de  la  zone  cible  est  illumine  et  les 
echos  de  tous  les  points  sont  recueillis  dans  la  mosaique  radar.  Le  traitement  de 
largeurs  de  faisceaux  finies  est  toutefois  congu  pour  ameliorer  le  realisme  physique  en 
tenant  compte  du  diagramme  d’antenne  reel.  Les  resultats  initiaux  montrent  que,  selon 
la  largeur  du  faisceau  efficace,  on  peut  reduire  le  nombre  de  boucles  de  traitement  de 
jusqu’a  50  %,  tout  en  maintenant  une  bonne  qualite  d’image  pour  ce  qui  est  de  la 
reponse  de  cible  reconstruite  (largeur  d’impulsion).  Les  ameliorations  possibles  pour 
le  traitement  RSO  3D  sont  encore  plus  importantes. 
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Executive  summary 


Typical  synthetic  aperture  radar  (SAR)  images  are  two-dimensional,  providing  range 
and  azimuth  information,  but  furnish  few  details  with  respect  to  elevation.  In  the  case 
of  Through- Wall  Sensing  (TWS)  radar,  having  a  third  imaging  dimension  could  be 
easily  implemented  and  be  of  immediate  use  in  target  detection  and  identification. 

This  work  outlines  one  approach  to  extend  SAR  to  three-dimensional  imaging:  replace 
the  single  antenna  element  by  an  array  of  antenna  elements.  The  simplest  such 
implementation  is  a  vertically  oriented  linear  array,  and  the  ensuing  transverse  motion 
of  the  assembly  produces  data  equivalent  to  that  of  an  extended  synthetic  planar  array. 
Image  reconstruction  proceeds  as  before,  where  range  resolution  is  achieved  via  pulse 
compression,  azimuth  resolution  via  angular  diversity  of  multiple  ’’look”  locations,  and 
the  new  dimension,  elevation,  derives  from  the  vertical  extent  of  the  physical  array. 
Rather  than  frequency-domain  based  reconstruction,  the  time-domain  method  of 
backprojection  was  selected  for  this  application.  It  is  simple  to  implement,  has 
localized  sidelobes,  and  it  can  handle  path  adjustments  for  motion  compensation  and 
propagation  through  a  dielectric  layer  such  as  a  building  wall. 

This  report  contains  two  sections,  the  first  outlining  data  and  image  processing  for  a 
three-dimensional  application,  and  the  second  part  details  finite  beamwidth  processing 
for  the  time-domain  backprojection  method.  Standard  backprojection  implicitly 
assumes  omnidirectional  antennas,  whereby  the  entire  target  area  is  illuminated  and 
returns  are  collected  from  all  points  in  the  imaging  grid.  This  decreases  the  efficiency 
of  the  algorithm  by  summing  over  unnecessary  grid  points.  Finite  beamwidth 
processing,  however,  is  designed  to  be  more  physically  realistic  by  accounting  for  the 
actual  antenna  pattern. 

This  work  shows  that  reducing  the  effective  beamwidth  for  2D  can  reduce  the  number 
of  processing  loops  by  up  to  50%,  while  still  maintaining  good  image  quality  in  terms 
of  the  reconstructed  target  response  (pulse  width).  Potential  improvements  in  3D  SAR 
processing  are  even  greater.  Simulations  carried  out  in  this  report  indicate  there  will  be 
an  optimum  beamwidth:  the  commonly  used  omnidirectional  configuration  provides 
good  results  at  the  cost  of  excessive  iterations  while  a  narrow  beamwidth  is  more 
efficient,  but  the  image  quality  is  impaired. 
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Sommaire 


Les  images  typiques  de  radars  a  synthese  d’ouverture  (SAR)  sont  bidimensionnelles  et 
foumissent  des  donnees  de  distance  et  d’azimut,  mais  peu  de  details  sur  le  site.  Dans  le 
cas  du  radar  de  visualisation  a  travers  les  murs,  une  configuration  d’imagerie  3D  serait 
facile  a  mettre  en  oeuvre  et  pourrait  etre  utilisee  immediatement  pour  la  detection  et 
P  identification  de  cibles. 

Les  presents  travaux  indiquent  comment  aj outer  une  troisieme  dimension  a  fimagerie 
SAR  en  remplagant  l’element  d’antenne  unique  par  un  reseau  d’elements  d’antennes. 
La  configuration  la  plus  simple  serait  un  reseau  lineaire  oriente  verticalement,  dont  le 
mouvement  transversal  produit  des  donnees  equivalant  a  celles  d’un  reseau  plan 
synthetique  etendu.  La  reconstruction  d’ image  se  fait  comme  auparavant;  la  resolution 
en  distance  etant  obtenue  par  compression  d’ impulsions  et  la  resolution  en  azimut  par 
la  diversite  angulaire  de  multiples  lieux  de  visee;  la  troisieme  dimension,  le  site,  est 
obtenue  a  partir  de  l’etendue  du  reseau  physique  suivant  la  verticale.  Pour  cette 
application,  la  methode  de  retroprojection  dans  le  domaine  temporel  a  ete  retenue  de 
preference  a  la  reconstruction  dans  le  domaine  des  frequences.  Le  reseau  est  simple  a 
installer  et  a  des  lobes  lateraux  localises;  il  permet  de  regler  les  trajets  de  maniere  a 
compenser  le  mouvement  et  gere  la  propagation  des  ondes  a  travers  une  couche 
dielectrique  telle  qu’un  mur  de  batiment. 

Le  rapport  se  divise  en  deux  parties;  la  premiere  porte  sur  le  traitement  des  donnees  et 
des  images  dans  une  application  tridimensionnelles,  et  la  seconde  examine  en  detail  le 
traitement  de  largeurs  de  faisceau  finies  selon  la  methode  de  retroprojection  dans  le 
domaine  temporel.  La  methode  de  retroprojection  courante  suppose  l’usage  d’antennes 
omnidirectives,  au  moyen  desquelles  l’ensemble  de  la  zone  cible  est  illumine  et  les 
echos  de  tous  les  points  sont  recueillis  dans  la  mosa'fque  radar.  Cette  methode  diminue 
Pefficacite  de  Palgorithme  en  faisant  la  sommation  de  points  inutiles.  Le  traitement  de 
largeurs  de  faisceaux  finies  est  toutefois  congu  pour  ameliorer  le  realisme  physique  en 
tenant  compte  du  diagramme  d’antenne  reel. 

Les  presents  travaux  montrent  qu’en  diminuant  la  largeur  efficace  du  faisceau  pour  le 
2D,  on  peut  reduire  le  nombre  de  boucles  de  traitement  de  jusqu’a  50  %,  tout  en 
maintenant  une  bonne  qualite  d’image  pour  ce  qui  est  de  la  reponse  de  cible 
reconstruite  (largeur  d’impulsion).  Les  ameliorations  possibles  pour  le  traitement  SAR 
3D  sont  encore  plus  importantes.  Les  simulations  effectuees  dans  le  cadre  des  travaux 
prevoient  une  largeur  de  faisceau  optimale  :  la  configuration  omnidirective 
couramment  utilisee  produit  de  bons  resultats  au  prix  d’ iterations  excessives,  alors 
qu’une  faible  largeur  de  faisceau  s’avere  plus  efficace,  la  qualite  de  l’image  etant 
toutefois  degradee. 


Barrie,  G.  2004.  Through-Wall  Synthetic  Aperture  Radar  (TWSAR)  3D  Imaging: 
Algorithm  Design.  DRDC  Ottawa  TM  2004-224.  R  &  D  pour  la  defense  Canada  - 
Ottawa. 
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1.  Introduction 


Typical  SAR  images  are  two-dimensional,  providing  range  and  azimuth  information, 
but  furnish  little  detail  with  respect  to  elevation.  The  underlying  principle  is  that 
transverse  motion  of  a  single  (point)  antenna  element  generates,  or  synthesizes ,  a  data 
set  equivalent  to  that  produced  by  an  extended  linear  array. 

In  three-dimensional  imaging,  we  can  readily  extend  the  concepts  from  two- 
dimensional  SAR  by  increasing  the  dimensionality  of  the  generating  element  by  one. 
That  is,  the  single  zero-dimensional  element  is  replaced  with  a  one-dimensional  linear 
array.  It  follows  that  transverse  motion  of  this  array,  oriented  perpendicular  to  the 
direction  of  motion,  produces  data  approximating  an  extended  two-dimensional  planar 
array.  The  ensuing  reconstruction  proceeds  as  before  where  range  resolution  is 
achieved  via  pulse  compression,  azimuth  resolution  via  angular  diversity  of  multiple 
’’look”  locations,  and  the  new  dimension,  elevation,  derives  from  the  vertical  extent  of 
the  physical  array. 

This  report  contains  two  sections,  the  first  outlining  data  and  image  processing  for  a 
three-dimensional  application,  and  the  second  part  details  finite  beamwidth  processing 
for  the  time-domain  backprojection  method.  Standard  backprojection  implicitly 
assumes  omnidirectional  antennas,  whereby  the  entire  target  area  is  illuminated  and 
returns  are  collected  from  all  points  in  the  imaging  grid.  This  decreases  the  efficiency 
of  the  algorithm  by  summing  over  unnecessary  grid  points.  Finite  beamwidth 
processing  is  designed  to  be  more  physically  realistic  by  accounting  for  the  actual 
antenna  pattern. 


Within  the  context  of  radar  imaging,  we  can  consider  three  different  coordinate 
systems,  defined  as  the  object,  data  and  image  space,  portrayed  in  Figure  1. 


a),  object  space:  g(x,y)  b).  data  space:  d(x,t)  c).  image  space:  f(x,y) 


Figure  1.  Coordinate  systems  used  in  synthetic  aperture  imaging. 

(Images  are  for  2D;  extension  to  3D  is  described  in  text.) 

The  object  space  g(x,y,z),  is  the  actual  physical  volume  (spatial  domain)  illuminated 
by  the  radar.  The  results  of  SAR  data  acquisition  and  pulse  compression  are  captured 
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in  the  data  space,  d(x,z,t) .  The  range  value  R  and  return  time  delay  t  are 
interchangeable  via  the  transformation  t  =  2R/  c.  In  moving  from  2D  to  3D,  the 
components  that  make  up  R  are  modified  to  include  elevation-dependent  range,  for 

V2  2  /  2  2  2 

x  +  j/  vs .R3D=^x  +  y  +  z  . 

In  data  space,  the  echoes  for  a  given  point  target  trace  out  a  hyperbola  due  to  the 
motion  of  the  radar.  The  two-dimensional  function  representing  this  “spreading”  of 
signal  energy  is  called  the  point  spread  response  (PSR).  Energy  contained  within  the 
PSR  must  be  collected,  or  focused  back  to  the  original  point  function,  within  the 
domain  of  image  space,  f(x,  y ,  z) . 

The  extension  from  two  to  three  dimensions  is  straightforward: 


2D 

3D 

g(x,y) 

g(x,y,z) 

d(x,t) 

d(x,  z,t ) 

f{x,y) 

f(x,y,z) 

2 
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2.  Three-dimensional  SAR 


2.1  System  geometry 

Here  x  denotes  cross-range,  or  azimuth ,  y,  the  range  coordinate  and  z  is  elevation  (refer 
to  Figure  2).  By  collecting  and  processing  radar  data,  one  attempts  to  reconstruct  an 
image  of  objects  contained  within  this  area.  For  a  simple  system  model,  we  consider  a 
discrete  set  of  point-reflectors  with  reflectivity  crp ,  (p  =  1,2,  •  •  •  ),  at  coordinates 

(wJ- 


Figure  2.  3D  imaging  geometry. 


2.2  Array  components 

The  linear  array  contains  M  elements,  each  located  at  position  vector  rm  where 
rm  —  m  =  1,2.  ..,M  .  The  w,v  and  w  are  cartesian  coordinates.  Further, 

there  are  a  number  of  synthetic  array  positions  n  =  1,2...,N  over  which  the  linear 
array  traverses,  resulting  in  Mx  N  antenna  element  coordinates1. 


1  It  is  difficult  to  find  a  clear  and  consistent  way  to  develop  the  notation  required  for  this  discussion,  but 
Soumekh's  text  [1]  provided  a  good  starting  point. 
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For  example,  an  M  =  5  element  linear  array  moving  across  TV  =  100  synthetic  array 
positions  would  generate  in  the  processor  memory  a  5  x  100  element  matrix. 
Consequently,  the  position  vector  of  the  mth  array  element  at  the  nth  synthetic  aperture 

position  is  rm  ={um,vm,wm). 


2.3  SAR  processor 

Backprojection  is  a  correlation  implementation  of  matched- filtering.  A  detailed 
description  of  the  backprojection  method,  along  with  implementation  details  can  be 
found  in  [2].  Presently,  it  suffices  that  the  idea  is  to  correlate  data  collected  at  each 
aperture  position  as  a  function  of  round-trip  delay  time.  This  means  the  correlation  is 
performed  with  respect  to  the  point  spread  response  (PSR)  obtained  in  the  space-time 
data  domain. 

Backprojection  sums  the  sampled  radar  returns  for  each  array  element  (pixel)  of  the 
image  map.  The  signal  obtained  at  each  aperture  position  is  time-shifted  to  match,  or 
align,  with  a  particular  pixel  element  in  the  image  map.  Following  this,  responses 
across  all  aperture  positions  are  combined  in  the  image  space. 


2.4  Channel  data 


For  3D  imaging,  the  target  echo  reaching  the  receiver  will  be  delayed  due  to  the  time 
of  propagation  between  the  Tx/Rx  array  and  the  target.  As  specified  earlier,  “targets” 
characterized  by  reflectivity  crp  are  located  at  positions  [xp 9yp9z  while  the  array 

element  is  located  at  (w,v,  w) .  Application  of  this  yields  the  following  signal  from 
scattering  point  1 : 


f 

st(r;t)  =  at  A  t 

V 


_M)2  +hi _v)2  +(zi  -w)2 


(i) 


which  is  a  function  of  antenna  element  position,  and  implicitly  is  a  function  of  time,  t 
as  a  consequence  of  propagation  delay.  In  this  expression,  A{t  -  r)is  the  impulse 
response  of  the  equipment  in  the  transmit  path,  including  additional  system 
impairments. 

The  composite  signal  received  at  a  given  array  element  is  the  summation  of  responses 
from  all  scattering  points  (see  Figure  3)  — 


4 1  ~~44P  M)2  +  - v)2  +  VP  ~ w)2 

p  v  c 

=  sx{r\t)  +  s2{r-,t)+... 


(2) 


Note:  To  keep  the  discussion  simple,  expressions  derived  here  are  for  the  monostatic 
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case.  To  consider  bistatic  configurations  as  depicted  in  the  diagrams,  replace  2R/c  by 
(RTx  +  Rrx)/c .  The  bistatic  case  was  implemented  in  computer  code. 


As  the  radar  array  traverses  the  synthetic  aperture,  signals  containing  combined  echoes 
from  all  targets  are  received  at  each  aperture  position.  The  collection  of  data  over  all 
element  positions  defines  the  data  space.  An  example  of  simulated  channel  data  is 
shown  in  Figure  4.  The  column  on  the  left-hand  side  shows  the  echos  received  at 
various  points  along  the  linear  array.  Each  antenna  is  referred  to  as  a  channel,  and 
there  are  nine  elements  in  this  particular  array.  The  received  data  from  odd  numbered 
channels  is  shown.  The  column  on  the  right-hand  side  shows  data  that  has  been 
collected  together  from  multiple  aperture  positions.  The  one-dimensional  data  set 
corresponding  to  frame  number  41  is  indicated  by  the  arrow.  Combining  multiple 
frames  results  in  the  two-dimensional  azimuth  vs.  range  delay  images  as  shown. 


DRDC  Ottawa  TM  2004-224 


5 


Figure  4.  Channel  data  for  scattering  centers  at  (2,2,0),  (5,5,2)  and  (8,8,4).  Units  are  in  meters. 


2.5  Image  formation 

Energy  contained  within  the  curves  must  be  collected,  or  focused  back  to  the  original 
point  function.  The  “optimum”  focusing  of  synthetic  aperture  data  is  obtained  by 
correlating  collected  data  s(f ;  t )  with  the  SAR  signature  at  a  given  grid  point, 

(.y. , y  . ,  ) .  In  effect,  we  require  a  space-time  domain  matched  filter  with  respect  to 

5(f ;  t ) ,  yielding  a  measure  of  the  reflectivity  at  that  point: 


2  / - 

fix, , y} , z k)  =  J | s(r; t)  A  t - -){x,  - uf  +  (v;  - v)2  +  (zk  - wf 


dt  dr.  (  3  ) 


where  r  is  the  position  vector  of  the  antenna  receive  elements  as  shown  in  Figure  2  on 
page  3.  Approximating  the  SAR  signature  as  a  delta  function,  — »  £(•),  we  obtain 

f{xi,yj,zk)  =  -uf  +(y,  -v)2  fzk  -  wf  j  dr 

=  \s(r;tyk{r))dr 
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The  time  convolution  is  calculated  as  — 


^^ijAr))=  s{r-,t)*  8(-t) 


=  J ' s(r;t) si t  - -  J(xt  - ; u  f  +  {y.  -  v)2  +  (zk  -  wf 

t  v  c 

=  s(V;  -  uf  +  (y,  -  v)2  +  (zt  -  wf 


dt 


(5) 


Returning  to  the  reflectivity  function,  the  full  image  is  obtained  by  calculating 
/U  ,  y j ,  zk  )  for  each  grid  point  in  the  image  space.  In  practice,  the  correlation 

integral  is  not  computed  directly.  Due  to  the  incremental  nature  of  both  aperture 
positions  and  time  samples,  a  two-dimensional  discrete  summation  is  performed. 
Further,  a  portion  of  the  calculation  requires  interpolation  since  the  recorded  time 
samples  do  not  correspond  exactly  to  grid  locations. 

Figure  5  is  a  representation  of  discrete  targets  producing  point  spread  responses 
(PSR’s)  in  the  space-time  domain.  This  diagram  shows  why  the  resolution  in  elevation 
is  so  limited.  In  comparison  with  horizontal  displacement,  there  are  only  a  limited 
number  of  elements  along  the  vertical  array,  and  a  correspondingly  reduced  number  of 
points  to  focus.  An  extended,  albeit  unrealistic,  array  10m  in  length  would  yield  a  PSR 
surface  as  in  Figure  6. 


Figure  5.  PSR’s  for  scattering  centers  at  (2,2,0),  (5,5,2)  and  (8,8,4).  Units  are  in  meters. 
Simulated  antenna  length  is  2m,  with  element  spacing  equal  to  0.25m. 
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Figure  6.  PSR  surface  for  extended  vertical  array. 


There  are  N  positions  in  the  synthetic  aperture,  and  M  positions  (channels)  in  the 
physical  array.  Grid  points  are  denoted  (xz- 9yJ9zk).  Therefore,  in  the  discrete 

computation,  we  will  have  a  total  of  five  nested  loops  to  perform  the  required 
summations  (I  x  J  x  K)  x  N  x  M. 


2.6  Positioning  errors 

In  an  ideal  case,  the  antenna  array  is  constrained  to  the  x-z  plane  and  therefore  has  no  y 
(range)  component  (that  is,  v  =  0).  In  general,  the  trajectory  followed  by  the  system 
will  not  be  as  precise  nor  will  it  be  as  predictable.  There  will  always  be  some  error  in 
the  presumed  positioning  of  the  Tx/Rx  array,  and  this  will  influence  the  precision  with 
which  images  can  be  reconstructed.  An  operational  system  will  have  methods  put  in 
place  to  track  and  compensate  for  these  deviations. 

In  terms  of  modeling  the  effects  of  positioning  errors,  these  could  be  accounted  for  by 
writing  the  generalized  position  vector: 

rs  —r  +5  =  (u  +  Su9  v  +  Sv9  w  +  dw) ,  (  6  ) 

and  defining  random  position  errors.  This  approach  to  predicting  system  performance 
has  been  carried  out  for  the  two-dimensional  case  [3],  but  can  be  generalized  to  the 
three-dimensional  case  by  recognizing  there  are  five  degrees  of  freedom  to  the  array's 
spatial  orientation:  the  position  vector  for  the  center  of  mass  and  the  orientation  angles 
for  azimuth  and  elevation  (Figure  7): 
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(w/>z/)  F  e  v 

<pt  0  <(pt<  2  n 
9i  0<0i<7r 

where  V  is  the  imaging  volume.  This  is  sufficient  to  characterize  the  array  since  each 
of  the  elements  are  constrained  to  lie  along  a  line  with  fixed  element  spacing.  For  the 
remainder  of  this  report  however,  we  assume  the  element  positions  are  known 
accurately. 


Figure  7.  Random  orientation  of  linear  array.  Center  of  mass  (CM)  corresponds  to  origin. 
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3.  Backprojection techniques 


3.1  Standard  backprojection 

A  major  concern  in  synthetic  aperture  imaging  application  of  backprojection 
techniques  involves  the  high  computational  effort.  For  each  element  of  the  synthetic 
aperture,  one  is  forced  to  consider  the  contribution  from  each  grid  point  within  the 
imaging  area  or  volume. 


To  perform  the  required  processing,  the  algorithm  requires  position  information  for 
both  Tx  and  Rx  elements  as  well  as  all  of  the  imaging  grid  points  (refer  to  Figure  8). 
This  allows  calculation  of  forward  and  reverse  paths,  resulting  in  a  value  for  the  2-way 
path,  which  is  ultimately  rescaled  in  terms  of  the  time  delay. 


(ij) 

-A 

IJ 

Figure  8.  Geometry  for  range  calculation  based  on  nth  array  element  and  grid  point  (i,j). 

The  backprojection  algorithm  then  matches  this  calculated  delay  against  recorded  data, 
either  voltage  or  field  strength  vs.  time,  interpolating  the  amplitude,  A ”  received  at 

that  grid  point  (Figure  9).  This  is  generally  an  interpolated  value  because  delay  times 
rarely  correspond  to  precise  sampling  points. 
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Figure  9.  Amplitude  interpolation  based  on  time  delay  for  sensor  position  n  and  grid  coordinates  (ij). 

To  clarify  notation,  the  66  A  ”  indicates  this  is  an  estimated  data  value  obtained  from  the 
nth  channel  (aperture  position),  and  is  allocated  to  the  grid  position  (i,j). 

A  straightforward  generalization  of  this  for  3D  imaging  is  A™  .  In  this  case,  n  remains 

the  synthetic  array  position  while  m  denotes  the  physical  array  element  position  along 
the  vertically  oriented  linear  array.  Grid  positions  for  the  three-dimensional  image 
space  are  labelled  (i,j,k). 

The  problem  with  this  technique  is  that  a  certain  amount  of  energy  is  spread,  or 
dispersed,  throughout  the  imaging  grid,  attributed  to  points  from  which  it  did  not 
originate.  As  a  consequence,  one  particular  field  value  is  applied  to  a  number  of 
points,  as  indicated  in  Figure  10.  Consider  a  situation  where  the  time  increment  is 

±  At .  The  value  A ”  is  then  applied  to  all  grid  points  falling  within  the  range  r..  ±  Ar , 

creating  artifacts  along  this  band.  This  non-uniqueness  is  a  characteristic  of  inverse 
problems  where  a  distinct  input  does  not  result  in  a  distinct  output-there  exist  a  range 
of  values  that  satisfy  the  initial  conditions. 
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Figure  10.  Schematic  of  image  grid  assignment  based  on  interpolated  data  (omnidirectional  antenna). 

This  is  mitigated  somewhat  by  the  principle  of  constructive  interference  —  the  grid 
points  corresponding  to  persistent  targets  will  continually  get  reinforced  with  each 
aperture  position,  increasing  the  SNR,  while  other  influences  fade  into  the  background 
noise. 

3.2  Finite  beamwidth  processing  for  2D 

Implicit  in  this  approach  is  the  assumption  that  the  antennas  are  isotropic.  In  a  strip 
map  application  such  as  this,  the  integration  angle  approaches  180°,  and  in  the  limit, 
approaches  an  infinite  aperture.  This  is  certainly  the  easiest  implementation  since 
loops  over  grid  coordinates  are  unrestricted.  A  more  realistic  situation  is  portrayed  in 
Figure  11,  where  the  antenna  possesses  a  finite  beamwidth,  0. 

In  this  situation,  the  shaded  area  represents  the  area  of  illumination,  where  the 
beamwidth  is  set  arbitrarily  at  say,  the  3  or  10  dB  points.  The  purpose  of  this  is  that 
points  outside  of  the  shaded  area  are  neither  illuminated,  nor  do  they  produce  a  radar 
return.  There  is  still  some  degree  of  non-uniqueness,  but  the  reduced  beamwidth 
reduces  its  extent. 
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Figure  11.  Finite  beamwidth  grid  point  assignment. 


To  implement  this,  we  require  a  method  to  limit  the  grid  points  to  those  that  lie  within 
the  illumination  area.  For  the  particular  case  considered  here,  namely  a  strip-map  SAR 
routine,  the  beam  axis  is  always  oriented  perpendicular  to  the  path  of  motion.  The 
more  general  case  as  depicted  in  Figure  1 1  is  more  appropriate  for  squint  or  spotlight- 
SAR  applications,  which  are  not  considered  in  this  report.  To  work  through  the  details 
of  the  algorithm,  refer  to  Figure  12. 


nominal 


First,  note  the  grid  spacing  for  azimuth  and  range  respectively,  Ax,  Ay  which  are  used 
as  the  basic  units  of  distance  in  the  following  text.  Next,  the  nth  synthetic  aperture 
element  is  located  at  [xn ,  yn  ) .  The  trajectory  of  the  radar  platform  can  be  displaced 
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from  the  nominal  path,  and  this  is  referred  to  as  the  “offset.”  Offset  is  denoted  simply 
as  Sr  =  yn  ,  and  this  value  will  be  positive  or  negative  depending  on  its  position 
relative  to  the  nominal  path.  By  definition,  the  range  of  the  imaging  grid  extends  from 
y  =  0  to  ymax.  The  offset  is  therefore  used  to  account  for  path  deviation  off  of  the  y- 
axis.  The  outer  loop  is  for  range  values,  moving  fromy  =  0  out  to  the  maximum  range. 
At  each  synthetic  aperture  (SA)  position,  we  look  downrange  to  grid  element  j  (that  is, 
along  the  perpendicular).  The  total  range  from  SA  element  to  grid  point  j  is  therefore 
Sr  +  jAy  (refer  to  Figure  12). 

We  have  specified  the  antenna  half-beamwidth  as  a  =  0/2  ,  and  can  now  calculate 
the  azimuth  extent  for  a  given  range,  LJX : 

Ljx  =  (Sr  +  y*Ay)tan<x  (  7  ) 

This  value  now  has  to  be  converted  into  upper  and  lower  index  values  for  the  grid 
points.  To  place  LJx  in  its  proper  location  in  the  imaging  grid,  we  need  the  x- 

component  of  the  nth  aperture  position,  xn .  To  obtain  the  lower  (or  left-hand)  limit, 

subtract  the  azimuth  segment  from  xn  and  divide  by  the  grid  spacing,  Ax  .  To  stay 
within  the  illuminated  area,  take  the  upper  integer  bound.  For  the  upper  (or  right-hand) 
limit,  we  add  the  segment  to  xn  and  take  the  lower  integer  bound.  The  loop  across  x 
grid  coordinates  at  a  given  range  increment  /,  for  synthetic  aperture  n  is  therefore: 

Ax 

Contrast  this  to  the  routine  without  modification.  In  that  situation,  the  loop  runs  from 
\<i<  DX . 

The  following  graphs  depict  a  progression  of  imaging  grids  corresponding  to  different 
locations  for  the  radar  platform  (Figure  13).  From  top  to  bottom  of  each  column,  the 
radar  moves  from  the  upper  to  lower  left-hand  comer.  There  are  five  aperture  positions 
in  total,  which  is  artificially  low,  but  still  serves  to  compare  focusing  of  targets-there 
are  two  in  the  image,  somewhere  along  each  arc.  At  this  time  in  the  reconstmction 
process,  the  azimuth  remains  undetermined.  Moving  from  left  to  right,  the  effective 
beamwidth,  and  therefore  the  illumination  area,  is  reduced  from  that  of  an 
omnidirectional  (-360°)  antenna  to  approximately  45°. 


*9+K 

Ax 


(8) 
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45° 


Omni 


120° 


Figure  13.  Image  maps  corresponding  to  discrete  radar  platform  positions. 


As  successive  images  corresponding  to  individual  element  positions  are  added,  the 
grid  points  where  arcs  overlap  tend  to  be  reinforced  while  other  regions  fade  (Figure 
14).  A  major  component  of  image  processing  therefore,  is  constructive  interference. 
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Omni  120° 


Range  [m]  Range  [m] 


Figure  14.  Reconstructed  images  corresponding  to  Figure  13. 

In  the  table  below,  loop  efficiency  reflects  the  reduction  in  the  number  of  iterations 
required  to  work  through  the  imaging  grid  as  in  Figure  14.  The  image  quality  criterion 
is  the  full- width  at  half-maximum  (FWHM). 


Table  1.  Comparison  of  processing  efficiency 
(Target  at  range  2.5m,  Synthetic  aperture  2.0m) 


Beamwidth 

Loop  efficiency  [%] 

FWHM  [m] 

omni 

0.0 

1.07 

120° 

12.8 

1.01 

CO 

o 

o 

22.3 

1.12 

cn 

o 

51.6 

1.25 

3.3  Restricted  beamwidth  processing  for  3D 

Having  established  this  method  for  the  2D  case,  we  can  now  (with  suitable 
generalizations)  apply  it  to  three-dimensional  imaging.  The  main  principle  remains  the 
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same:  take  one  step  downrange  and  adjust  the  grid  point  indices  to  account  for  the 
antenna  beam  pattern.  For  the  3D  case,  instead  of  determining  segment  length  along 

the  azimuth,  Ljx ,  we  now  consider  two  segments  in  the  x-z  plane,  namely  Ljx  and  Ljz . 
Once  again,  j  is  the  index  for  the  range  coordinate. 

To  account  for  the  volume  illuminated  by  the  antenna,  we  specify  two  angles  where 
0  is  the  horizontal  beamwidth  and  ®  is  the  vertical  beamwidth.  This  leads  to  the 
half-beamwidth  angles  a  =  0  /  2  and  /?  =  ®  /  2 .  In  turn,  this  results  in  — 


L'x  =  [Sr  +  /Ay  )  tan  a 

(9) 

L[  =  (Sr  +  y'Ay)tan  [3 

(10) 

These  quantities  are  shown  in  Figure  15  as  projections  on  the  xz  plane  at  a  given  range. 
In  general  the  shaded  area  will  be  elliptical,  but  to  simplify  calculations  the  projection 
is  assumed  to  be  rectangular. 

z 


Figure  15.  Projection  of  beam  in  3D  coordinate  system. 


Following  the  same  reasoning  as  before,  the  indices  for  x-  andz-  grid  points  are 
restricted  to: 
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(11) 


Ax 


*"+4 

Ax 


and 


zm+L{ 

•  k  • 

1 

N 

3 

1 

_ 1 

Az 

.  rv 

Az 

To  consider  a  more  accurate  depiction  of  how  the  antenna  beam  projects  onto  the  xz 
plane,  one  could  consider  an  ellipse  (Figure  16),  as  it  would  account  for  both 
horizontal  and  vertical  beamwidths. 

z 


X 


Figure  16.  Interpolation  of  grid  points  for  elliptic  projection. 

However,  it  would  be  necessary  to  calculate  grid  points  as  an  implicit  function  of  p , 
and  Az.  The  problem  stems  from  the  fact  that  p  will  change  depending  on  the  value  of 
z  =  nAz.  To  make  matters  worse,  you  need  to  know  the  extent  of  the  major  and  minor 
axes  to  obtain 


P  = 


b 2 

l-e2  cos2  y 


(13) 


where  e 2  =  1  -b2 ! a2 . 
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To  determine  the  horizontal  extent  of  the  ellipse,  and  consequently  the  min/max  loop 
indices,  we  would  have  to  determine 


L  = 


l-e  cos  y J 


-  ( nAz )2 


(14) 


for  each  elevation  increment.  This  would  cause  excessive  cpu  time  and  it  is  suggested 
that  the  xz  projection  remain  rectangular. 
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4.  Conclusion 


In  this  report,  two-dimensional  synthetic  aperture  radar  (SAR)  is  generalized  for  three- 
dimensional  imaging  application,  for  both  data  and  image  processing  operations.  The 
key  to  effective  algorithm  implementation  is  having  a  clear,  consistent  description  of 
the  required  processing  steps — these  have  been  illustrated  graphically  and  with  the 
governing  equations. 

Rather  than  frequency-domain  based  reconstruction,  the  time-domain  method  of 
backprojection  was  selected  for  this  application.  This  is  simple  to  implement, 
generates  localized  sidelobes,  and  it  can  handle  path  adjustments  for  motion 
compensation  and  propagation  through  a  dielectric  layer  such  as  a  building  wall. 

Another  advantage  of  backprojection  is  that  intermediate  stages  of  computation  are 
intuitive  and  easily  visualized  since  there  are  no  time-frequency  transformations.  For 
example,  resolution  in  elevation  is  limited  by  the  physical  length  of  the  linear  array, 
but  another  view  of  how  these  limitations  arise  was  presented  based  on  the  focusing  of 
returns  in  3D  data  space,  d(x,  z,  t). 

As  this  is  currently  envisioned  as  a  ground-based  application,  positioning  errors  arising 
from  the  arbitrary  orientation  of  a  linear  array  moving  across  irregular  terrain  need  to 
be  considered.  A  full  analysis  was  not  undertaken,  but  the  basic  problem  was  outlined, 
along  with  a  description  of  how  the  situation  could  be  simulated.  This  will  take  on 
greater  significance  when  motion  compensation  is  considered  in  the  future. 

A  major  drawback  with  backprojection  is  the  severe  computational  requirement  of 
having  to  loop  over  every  grid  point  in  image  space.  A  method  to  restrict  the  area  or 
volume  by  adjusting  loop  indices  was  described.  This  results  in  a  smaller  area  of 
illumination  within  the  target  area,  producing  a  more  realistic  scenario.  Initial  results 
show  that  depending  on  the  width  of  the  effective  beam,  the  number  of  loops  can  be 
reduced  by  up  to  50%,  while  still  maintaining  good  image  quality  in  terms  of  the 
reconstructed  target  response  FWHM.  Potential  improvements  in  3D  SAR  processing 
are  even  greater. 

Simulations  indicate  there  will  be  an  optimum  value  for  the  beamwidth.  The 
omnidirectional  beam  commonly  used  provides  good  results  at  the  cost  of  excessive 
iterations.  A  narrow  beamwidth  is  more  efficient,  but  the  image  quality  is  impaired-a 
major  aspect  of  SAR  processing  is  building  up  an  image  by  viewing  it  from  many 
positions.  This  angular  diversity  is  enhanced  by  using  a  broad-beam  antenna,  so  the 
simulations  support  the  intuition  that  the  best  setting  will  be  somewhere  between  the 
two  extremes,  and  that  incorporating  the  antenna  pattern  into  the  backprojection 
method  will  result  in  more  effective  and  efficient  image  reconstruction. 
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